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ABSTRACT 

We investigate the relationship between 8.44 GHz brightness temperatures and 1.4 to 8.44 GHz radio 
spectral indices with 6.2 /zm polycyclic aromatic hydrocarbon (PAH) emission and 9.7 /zm silicate 
absorption features for a sample of 36 local luminous and ultra-luminous infrared galaxies. We find 
that galaxies having small 6.2 /im PAH equivalent widths (EQWs), which signal the presence of weak 
PAH emission and/or an excess of very hot dust, also have flat spectral indices. The three active 
galactic nuclei (AGN) identified through their excessively large 8.44 GHz brightness temperatures 
are also identified as AGN via their small 6.2 /xm PAH EQWs. We also find that the flattening of 
the radio spectrum increases with increasing silicate optical depth, 8.44 GHz brightness temperature, 
and decreasing size of the radio source even after removing potential AGN, supporting the idea that 
compact starbursts show spectral flattening as the result of increased free-free absorption. These 
correlations additionally suggest that the dust obscuration in these galaxies must largely be coming 
from the vicinity of the compact starburst itself, and is not distributed throughout the (foreground) 
disk of the galaxy. Finally, we investigate the location of these infrared-bright systems relative to the 
main sequence (star formation rate vs. stellar mass) of star-forming galaxies in the local universe. 
We find that the radio spectral indices of galaxies flattens with increasing distance above the main 
sequence, or in other words, with increasing speciflc star formation rate. This indicates that galaxies 
located above the main sequence, having high specific star formation rates, are typically compact 
starbursts hosting deeply embedded star formation that becomes more optically thick in the radio 
and infrared with increased distance above the main sequence. 

Subject headings: galaxies: active - galaxies:starbursts - infrared:galaxies - radio continuum:galaxies - 
stars:formation 



1. INTRODUCTION 

Distinguishing the physical processes driving the en- 
ergetics of infrared-bright galaxies is a non-trivial task. 
This is especially true for luminous infrared galaxies 
(LIRGs) whose infrared (IR; 8 — 1000 /im) luminosities 
exceed iiR > 10^^ Lq. Immense column densities of dust 
towards the nuclei in these galaxies typically render tra- 
ditional optical diagnostics (e.g. emission line flux ra- 
tios) difficult to interpret. In the redshift range span- 
ning 1 ^ z < 3, such infrared- luminous galaxies appear 
to be much more common and dominate the star forma- 
tion rate density, being an or der of magnitude larger than 
today (e.g., p iarv & Elba d [20?)ll: ILe Floc'h et all [20051: 
Caputi et alTIToOTI: iMagnelli et al.l 120111 : iMurphv et al.l 



201 lal) . 



At these epochs, as well as today, there seems to 
be a distinction between the modes of star formation 
in infrared-bright sources that emerges by comparing 
their star formation rates and stellar masses. The ma- 
j ority of sources lie along a seemingly redshift depen- 
dent (jNocskc ct al." 2007; Elbaz ct al. 2007; Daddi et all 
I2007t iPannclla ct al. 2009; Magdis et al. 201^) "main- 
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sequence" where galaxies appear to be undergoing a 
mode of distributed star formation. A minority of 
sources, having compact star formation exhibited by ele- 
vated star formation rate surface densities, are instead 
found to lie systematically above the main sequence, 
having a high star formation rate per unit stellar mass 
(lElbaz et al.l |2011l ) . The difference between these two 
populations is most likely related to the gas content 
and fueling of star formation in these systems, as illus- 
trated by a similar sep aration of the same sources in th e 
Schmidt law diagram (|Schmidtl [19591 : IKennicuttl[T998h . 
Main sequence galaxies follow the classical Schmidt law 
while compact starbursts have star formation rates that 
are an order of magnitude larg e r for the sarne gas sur- 
face density (|Daddi et al.|[20Ta [Gmiel et al.|[2010l ). If 
the modes of star formation that are responsible for plac- 
ing galaxies on and off the star-forming main sequence 
are the same at all redshifts, we should be able to use 
samples of local galaxies, for which we can obtain much 
more detailed information on small physical scales, to 
shed light on the detailed astrophysics of such distant 
systems. 

In the local universe, it is well-known that LIRGs 
and uhraluminous LIRGs (ULIRGs; Lir > lO^^i©) 
appear to be undergoing an intense starburst phase. 
Within these systems are compact star-forming re- 
gions that have been been triggered predominantly 



through major mer gers (see e.g.JArmus et al.lll987 
T989, 1990; Sandcr setaIlll988aU bl: IMurphv et al 
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pact cores of ULIRGs have been resolved by high (i.e., 
sub-arcsecond) resolution, mid - infrared (i.e., 8 — 25 ^m) 
imaging from the ground fe.g.. iSoifer et al.ll2000D . How- 
ever, given the large beams of typical space-based far- 
infrared (i.e., A > 25 /im) telescopes, which are able to 
directly measure the bulk of the re-radiated energy from 
compact starbursts and active galactic nuclei (AGN) 
without suffering significantly from extinction, interfero- 
metric radio/mm/submm observations are currently the 
best means to directly resolve the sizes of these highly 
energetic, and often extremely compact, sources (e.g., 
iLonsdale et"an '2006': 'Saka moto et al.ll2009[ ). 

iCondon et al., (1991) investigated the radio properties 
for a sample of 40 infrared-luminous galaxie s included in 
the IRAS revised Bright Galaxies Sample (jSoifer et al.l 
[19891 : ISanders et all [2001 . Each source has a 60 //m 
flux density larger than 5.24 Jy and a dust tempera- 
ture in the range of 60 < Td < 80 K, suggesting that 
their angular size at 60 /im must be >0'.'25. Thus, ra- 
dio maps made at such a resolution should res olve the 
far- infrared source, leading [Condon et al.l ()1991[ ) to map 
the sample with 0'.'25 resolution at 8.44 GHz using the 
Very Large Arr ay (VLA) . In this p aper we compare 
the results from ICondon et all (|1991[ ). which measured 
8.44 GHz sizes and brightness temperatures of these com- 
pact, infrared-bright galaxies, with mid-infrared spectro- 
scopic diagnostics taken as part of t he Great Observato - 
ries All-Sky LIRG Survey (GOALS: [Almus et al.l[2009[ ). 
In doing this we attempt to see how the mid-infrared 
spectral properties compare with the physical conditions 
of the starbursts (e.g., sizes and brightness temperatures) 
as measured by the radio data. 

The paper is organized as follows: In §2 we describe the 
existing radio and mid-infrared data used in the analysis. 
Then, in §3, we compare various mid-infrared spectral 
and radio continuum properties of the sample, which are 
then discussed in §4. Finally, we summarize our main 
conclusions in §5. 

2. DATA AND ANALYSIS 

The galaxies in this sample are drawn from the 
IRAS revised Bright Galaxies Sample (Soifer et al. 1989; 
ISanders et all [2003) having 60 /im flux densities larger 
than 5.24 Jy and far-infrared (FIR; 42.5 - 122.5/im) lu- 
minosities > 10^^'^^ Lp, . The 40 syste r ns tha t meet these 
criteria were imaged bv lCondon et al.l (|1991[ ) at 8.44 GHz 
with 0'.'25 resolution, sometimes resolving multiple com- 
ponents. Of the resolved galaxies, radio spectral indices 
measure d between 1.49 and 8 .44 GHz are available for 36 
sources (jCondon et al.l[T99ll ). which defines the sample 
used in the present analysis. We additionally make use 
of the derived 8.44 GHz sizes and brightness tempera- 
tures available for 28 of the sou rces, which are all given 
in Table [H As shown by Condon et al.l (|199H ). sources 
having 8.44 GHz brightness temperatures > lO'^'^K can- 
not be sustained by star formation alone, and must be 
powered by nuclear "monsters" (AGN). 

The mid-infrared spectral properties w ere collected 
by th e Spitzer Infrared Spectr ograph (IRS: iHouck et al.l 
I200I as part o f GOALS (lArmus et al I 
are taken from iStierwalt et al.l (|2013bl l; 
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sources, IRS observations were either not taken (MCG - 
03-12-002) or the slit missed the nucleus of the 
galaxy (IRAS F01417-M651, IRASF03359-I-1523, and 



IRASF17132-I-5313). The total IR (8 - 1000 /im) lu- 
minosities are taken from lArmus et al.l (|2009[ ). How- 
ever, for sources having multiple components as re- 
solved by the IRS, the fractional IR luminosity from 
each component is determined by the fraction of re- 
solved 24 /im emission. We include 6.2/im PAH equiv- 
alent widths (EQWs), 9.7Atm silicate optical de pths, and 
silicate strengths (see IStierwalt et al.l l2013bl [a[) in our 
analysis. The silicate strength at 9.7 /xm is defined as 
S9.7/im = log(/9.7/im/C9.7/im), where /9.7/im is the mea- 
sured flux density at the central wavelength of the ab- 
sorption feature and Cg.y^m is the expected continuum 
level in the absence of the absorption feature. While the 
0'.'25 radio maps are at a resolution that is significantly 
higher than the mid-infrared IRS data (i.e., 3'.'6), this 
should not affect our analysis or conclusions given that 
the central sources seem to dominate the luminosities of 
each system. Furthermore, the radio data themselves 
should be sensitive to extended emission on scales up to 
Rs5", which is larger than the re solution of the IRS data. 

Stellar masses are taken from lU et all (|2012D . who de- 
rived estimates by fitting optical through near-infrared 
spectral energy distributions, as well as assuming a 
fixed mass-to-li ght ratio and ff -band fiux den sities, for 
both Salpeter (jSalpeteii I1955D and Chabrier (jChabrieii 
[2OOI IMFs. The masses adopted here are the mean 
of these estimates aft er first converting them to a 
Kroupa (Kroupa 20011 ) IMF using the follow i ng sta - 
tistical relations as given by 'Bolzoncll a et al.l (|2010f l: 
logM^Chab) « logM4Salp)-0.23 and logM*(Krou) « 
logM*(Chab) -I- 0.04. For sources having multiple com- 
ponents, we di vided the tota l mass es up by the propor- 
tions given bv iHowell et al.l (|2010D . These authors de- 
rived stellar masses using Two Micron All Sky Survey 
(2MASS) Ks-hand (where possible) and IRAC 3.6 /im 
photometry, along with the z y mass-to-light ratios 
presented by iLacev et al.l (I2008D. We pre fer to use the 
mass estimates reported bv lU et all (|2012D as they have 
been more robustly determined. Each of these proper- 
ties is given in Table [1] Infrared luminosities are con- 
ve rted into a star form ation rate using the calibration 
of iMurphv erall (|2011bl ). which assumes a Kroupa IMF, 
such that, 
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We note that the coefficient here is 1.23 times larger, 
but within the scatter, of the empirically derived rela- 
tion between IR and radio free-free star formation rates 
(|Murphv et al.llMl) . 

3. RESULTS 

In this section we present the main findings by compar- 
ing the radio and mid-infrared properties of the sample. 

3.1. Radio and Mid-Infrared Spectral Properties 

In the top panel of Figure[T]we plot the 8.44 GHz linear 
sizes of each source against their radio spectral index, 
typically measured between 1.49 and 8.44 GHz (see Table 
[1]). While there is a large range in size among galaxies 
having radio spectral indices steeper than Q;f |gQ§| > 
0.5, sources with a^ ^gQ^^ ^ 0-5 all have small radio 
sizes. The mean and standard deviation in linear sizes 
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Table 1 







Radio 


and Infrared 


Properties 


of the Sample 


Galaxies 








Galaxy 


Dist.'' 


ilR/lO"'' 


^8.44 GHzb 
1.49 GHz 




Tb/lO**^ 6 


.2^lm EQW^ 






M./IO^O'* 




(Mpc) 


(^-0) 


(kpc2) 


(K) 


(/im) 






(A/©) 


NGC 34 


84.1 


3.09 


0.78 


33.96 


0.13 


0.45 


2.39 


-0.79 


2.08 


CGCG 436-030 


134.0 


4.90 


0.70 


89.38 


0.08 


0.35 


3.00 


-1.10 


1.20 


IRAS F01364-1042 


210.0 


7.08 


0.42 


25.41 


0.63 


0.39 


5.14 


-1.27 


2.96 


IRAS F01417-f 1651 


119.0 


4.37 


0.40 


15.64 


0.79 








1.17 


NGC 695 


139.0 


4.79 


0.88 " 






0.65 


0.57 


0.31 


9.77 


UGC 2369 


136.0 


4.68 


0.67 


12.76 


0.32 


0.57 


0.89 


-0.11 




IRAS F03359-I-1523 


152.0 


3.55 


0.31 


10.48 


0.63 










MCG -03-12-002 


138.5 


3.24 


0.65 


8.01 


0.63 










NGC 1614 


67.8 


4.47 


1.05 " 






0.61 


1.31 


-0.41 


2.47 


IRAS F05189-2524 


187.0 


14.45 


0.52 


43.51 


0.40 


0.03 


1.39 


-0.29 


9.80 


NGC 2623 


84.1 


4.37 


0.58 


25.03 


0.40 


0.27 


3.15 


-1.12 


1.71 


IRAS F08572-f3915 


264.0 


14.45 


0.27 


15.95 


0.79 


0.03 


6.96 


-3.58 


1.59 


UGC 04881 E 


178.0 


3.63 


0.69 


38.98 


0.25 


0.40 


2.00 


-0.81 


4.10 


UGC 05101 


177.0 


10.23 


0.59 


18.82 


3.16 


0.13 


3.65 


-0.78 


7.94 


IRASF10173+0828 


224.0 


7.24 


0.28 


< 7.63 


> 1.58 


0.35 


4.24 


-1.20 


2.55 


IRAS F10565-f 2448 


197.0 


12.02 


0.68 






0.51 


2.13 


-0.75 


3.49 


MCG -f 07-23-019 


158.0 


4.17 


0.55 


22.80 


0.40 


0.64 


2.31 


-0.55 




IRAS F1123H- 1456 


157.0 


4.37 


0.82 






0.60 


1.16 


-0.22 


7.04 


NGC 3690 E 


50.8 


3.63 


0.67 


11.55 


0.79 


0.38 


3.82 


-1.65 


0.98 


IRAS F12112-I-0305 


340.0 


22.91 


0.47 


50.27 


0.79 


0.30 


3.42 


-1.24 


5.87 


UGC 08058 


192.0 


37.15 


-0.06 


< 5.51 


> 79.43 


0.01 


> 1.84 


-0.48 


30.63 


VV 250a E 


142.0 


5.13 


0.64 


125.76 


0.08 


0.63 


1.80 


-0.67 


0.90 


UGC 08387 


110.0 


5.37 


0.64 






0.62 


3.01 


-1.01 


2.30 


NGC 5256 S 


129.0 


3.09 


0.64 " 






0.44 


1.50 


-0.47 


5.33 


UGC 08696 


173.0 


16.22 


0.63 


75.02 


0.63 


0.12 


4.24 


-1.37 


5.33 


IRASF14348-1447 


387.0 


24.55 


0.71 


149.04 


0.25 


0.25 


4.31 


-1.36 


6.79 


VV 340a N 


157.0 


4.37 


0.60 " 






0.58 


2.17 


-0.63 


4.61 


IRASF15163-I-4255 


183.0 


8.32 


0.78 


36.36 


0.25 


0.75 


1.06 


0.32 




IRAS F15250-I- 3608 


254.0 


12.02 


0.11 


< 7.51 


> 3.98 


0.03 


6.21 


-2.69 


2.16 


UGC 09913 


87.9 


19.05 


0.41 


10.09 


2.51 


0.17 


6.10 


-2.26 


4.48 


NGC 6286 S 


85.7 


1.62 


0.84 ° 






0.59 


2.28 


-0.40 


2.61 


IRAS F17132-f 5313 


232.0 


9.12 


0.67 


46.64 


0.13 








11.14 


IRAS F22491-1808 


351.0 


15.85 


0.41 


< 19.74 


> 0.79 


0.48 


3.23 


-1.04 


3.29 


NGC 7469 


70.8 


4.47 


0.68 


2.71 


1.26 


0.23 


0.70 


0.06 




IC 5298 


119.0 


3.98 


0.67 


5.78 


0.50 


0.12 


1.55 


-0.37 


4.15 


MRK 331 


79.3 


3.16 


0.66 


2.99 


0.79 


0.63 


1.59 


-0.35 


0.97 



^ Distances and total IR (8 — 1000 /im) luminosities arc taken from fArmus et all i200gr ■ 
^ Radio data are taken from [Condon et al.l ] 19911) . 

^ 6.2 fim PAH EQWs and silicate strengths (so.ry^m) arc taken from JStierwalt et aril2013lJl while silicate optical depths (ro.Ty^m) are take 
from lStierw alt et aj] J2013al) . 

^ The mean mass estimate given in lU et all i2012f) after converting to a Kroupa IMF. 

^ Spectal indices are based on measurements at 1.49 GHz and a frequency other than 8.44 GHz fsee lCondon et al.lll99lT) . 



among sources steeper and flatter than a^ ^gQiiz ~ 0.5 
are 0.12 ± 0.04 kpc and 0.18 ± O.lOkpc, respectively. 

By instead plotting the 8.44 GHz brightness tempera- 
ture against the radio spectral index (see bottom panel 
of Figure [Bi the r e is a general trend as pointed out 
by iCondon et al.l (|1991[ ) in which more compact (i.e., 
higher surface brightness) sources tend to have flatter 
radio spectral indices, consistent with compact sources 
becoming optically thick at low (i.e., i' < 5 GHz) radio 
frequencies. A similar trend is also seen when the ra- 
dio spectral index is compared to the mid-infrared size 
of the starburst as meas ured from the Spitzer IRS data 
(jDfaz-Santos et aTll2010t ). Galaxies with flatter spectral 
indices tend to have smaller fractions of extended mid- 
infrared (13.2 /im) continuum emission, again suggesting 
that the radio spectral index is related to the compact- 
ness of the starburst. 

In Figure [2] we plot the optical depth of the 9.7 /xm sil- 
icate feature against the radio spectral index. We find a 
correlation such that the silicate optical depth increases 
as the radio spectrum flattens. The only outlier in this 
trend is UGC 08058 (Mrk231), which has a very small 
silicate optical depth given its extremely flat (inverted) 



radio spectrum. The shallow silicate optical depth in 
UGC 08058 is most likely due to hot dust emission that 
fills in the absorption trough, thus the measured optical 
depth probably provides a lower limit. Additionally, its 
position in this diagram may also be complicated by the 
fact that the variable radio source in UGC 08058 is small 
enough in size (i.e., < Ipc) that its flat spectr um can be 
explained by synchrotron self- absorption (Co ndon et al.l 
[1991). If we instead compare the radio spectral indices 
with the silicate strength, which does not rely on a model 
fit and measures the apparent silicate optical depth di- 
rectly from the spectrum, we find a similar trend and 
UGC 08058 remains an outlier. 

3.2. Identifying Potential AGN-Dominated Sources 

In Figure [3] we look at the 6.2 /im PAH EQW, which 
measures the relative strength of the PAH flux with re- 
spect to the amount of hot dust emission at 6.2 /im, ver- 
sus radio spectral index. The 6.2 /im PAH EQW has 
been used to indicate the presence of an A GN since they 
typica.lly have very small PAH EQWs fe.g.. lGenzel et al.] 
Il998t lArmus et al.l [20071 ) . Specifically, starburst domi- 
nated systems appear to have 6.2 /xm PAH EQWs that 
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Figure 1. Top: The linear 8.44 GHz size of each source plotted 
against the radio spectral index (generally) measured between 1.49 
and 8.44 GHz (see Table [TJ. Bottom: The same as the top panel, 
except that the 8.44 GHz brightness temperature is now plotted 
against the radio spectral spectral index. As already shown by 
[Smdon et al. (199j), there is a general trend for sources having 
fiat spectral indices to be small in size and have large 8.44 GHz 
brightness temperatures. In both panels, sources with 8.44 GHz 
brightness temperatures larger than 10*-^ K, indicating the pres- 
ence of an AGN, are identified by filled circles. We additionally 
identify those sources categorized as being AGN dominated by 
their low 6.2 /xm PAH EQW (circles with crosses); note, each of 
the sources with an 8.44 GHz brightness temperature larger than 
lO^-S K also has a 6.2 /xm PAH EQW < 0.27 ^tm . 

are > 0.54 /.tm (|Brandl et all I2OO60 . while AGN have 
6.2 /xm PAH EQWs are < 0.27 /xm. We find that the 
galaxies having steep spectral indices (i.e., af;4gQ§| > 
0.6) never have extremely low 6.2//m PAH EQWs (i.e., 
<0.1/im). Additionally, all 3 sources identified as AGN 
via their excessively large 8.44 GHz brightness tempera- 
ture sit in the PAH-defincd AGN region. 

Using the 6.2 fim PAH EQW to split the sources into 
starburst and AGN dominated systems, we find that the 
mean radio spectral index for starburst dominated sys- 
tems is 0.74, with a standard deviation of 0.15. For the 
AGN dominated systems, we find a significantly fiatter 
mean spectral index of 0.45, albeit with a much larger 
scatter of 0.27. 

4. DISCUSSION 

Taking archival high resolution 8.44 GHz data from the 
literature, we have compared the radio and mid-infrared 
spectral properties for a sample of local (U)LIRGs. We 
now use these results to investigate how the dust is dis- 
tributed in these powerful, infrared-bright galaxies. We 
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Figure 2. The optical depth of the 9.7/.tm silicate feature plot- 
ted against the radio spectral index (generally) measured between 
1.49 and 8.44 GHz (see Table[TJ. Sources with 8.44 GHz brightness 
temperatures larger than 10^'^ K, indicating the presence of an 
AGN, are identified by filled circles. We additionally identify those 
sources categorized as being AGN dominated by their low 6.2 fim 
PAH EQW (circles with crosses); note, each of the sources with 
an 8.44 GHz brightness temperature larger than lO'*'^ K also has a 
6.2 fim PAH EQW < 0.27 ^tm. The dotted Une is an ordinary least 
squares fit to all galaxies while the dashed line excludes sources 
identified as possible AGN. Both fits illustrate a general trend for 
sources to have decreasing silicate optical depths as their radio 
spectra steepen. The source to the lower left, having a negative 
radio spectral index, is UGC 08058 (Mrk231) and is not included 
when fitting the data given that we only have a lower limit for its 
silicate optical depth. 
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-0.2 0.0 0.2 0.4 0.6 O.E 
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Figure 3. The 6.2 fira PAH EQW plotted against the radio spec- 
tral index (generally) measured between 1.49 and 8.44 GHz (see 
Table [T} . The horizontal lines indicate regions for which sources 
are thought to be AGN or star formation driven based on their mea- 
sured 6.2 fim PAH EQW. Sources having spectral indices steeper 
than Qi'lgQH^ 0.6 never appear to have extremely low 6.2/im 
PAH EQWs (i.e., <0.1 fim). Sources with 8.44 GHz brightness tem- 
peratures larger than 10'*'^ K, indicating the presence of an AGN, 
are identified by filled circles. The source to the lower left, having 
a negative radio spectral index, is UGC 08058 (Mrk231). 

additionally look at how the radio properties of these 
sources relate to their location in the star formation rate 
- stellar mass plane of star-forming galaxies. 

4.1. The Dust Distribution in Compact Starbursts 

In Figure [T] we have shown that the 8.44 GHz radio 
sizes (top panel) and brightness temperatures (bottom 
panel) of sources decrease and increase, respectively, as 
a galaxy's radio spectrum flattens. While it may not 
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be too surprising to find compact, flat-spectrum ACN, 
which indeed make up the two sources having the flattest 
radio spectra, we find that these trends persist even after 
identifying potential ACN-dominated sources, as indi- 
cated by low 6.2 /im PAH EQW (i.e., < 0.27 ^m) or high 
8.44 GHz brightness temperature [i.e., rb(8.44CHz) > 
lO^-^K]. 

Additionally, in Figure [21 we find a trend such that the 
silicate optical depth increases with the fiattening of the 
radio spectral index. Sources potentially dominated by 
ACN, as indicated by their low 6.2 fiin PAH EQW (i.e., 
< 0.27 ^m) or high 8.44 GHz brightness temperature [i.e., 
rb(8.44GHz) > lO'^-^K] are identified. An ordinary least 
squares fit to the data, excluding UGC 08058 for which 
we only have a lower limit on its silicate optical depth, is 
shown as a dotted line and yields the following relation; 



Tg.T^m = (-6.5 ± 1.0)a + (6.8 ± 0.7). 



(2) 



The fact that we see a trend with fiat spectrum sources 
having larger silicate optical depths lends support to the 
idea that starbursts show spectral fiattening as the re- 
sult of increased free-free absorption arising from more 
deeply embedded star formation. In fact, if we instead 
limit the fit to non-AGN sources, as indicated by ei- 
ther a low 6.2 fim PAH EQW or a high 8.44 GHz bright- 
ness temperature, the trend actually tightens; the resid- 
ual dispersion about the fit (see the dashed line in Fig- 
ure [2]), whose coefficients are consistent within errors 
[T9.7/im = (-4.9±1.0)a-h(5.6±0.6)],is 30% smaller than 
when fitting the entire sample. The increased dispersion 
by including AGN may not be surprising given that, un- 
like optically-thick starbursts, flat-spectrum AGN need 
not be heavily obscured. 

Since the radio observations are not significantly af- 
fected by dust extinction, we can assume that the ra- 
dio data are able to penetrate the starburst core itself. 
iCondon et al.l (|1991[ ) showed that the spectral fiatten- 
ing observed in the radio data can be explained by the 
starburst becoming increasingly compact and energetic. 
Thus, the correlation between the spectral flattening and 
increased silicate optical depth (Figure [2]) and the cor- 
relations between the spectral fiattening and the radio 
extent and brightness temperature (Figure [IJ most likely 
indicates that the dust obscuration in these sources oc- 
curs in the vicinity of the compact starburst itself, and 
not by extended dust located in the foreground galaxy 
disks. 

4.2. Radio Spectral Indices and the Main Sequence of 
Star-Forming Galaxies 

In Figure U] we plot the -g m ain sequence re- 
lation as given by lElbaz et al.l (|2Qllf) , which is estab- 
lished using a local sample of galaxies dete cted at Q^tm 
with the AKARI/Infrared Camera (IRC; lOnaka et al.l 
l2007f ). Most galaxies lie along a linear relation between 
star formation rate and stellar mass, such that the spe- 
cific star formation rate (sSFR — SFR/M*) is constant, 
being sSFR ~ 0.25 Gyr^^ with a dispersion (standard 
deviation) that is a factor of 1.82. We n ote that the 
z ^ main-sequence relation reported by lElbaz et al.l 
(|2007D has a slope that is slightly flatter than unity, be- 
ing 0.77. If we instead adopt this relation for the z ^ 
main sequence, our conclusions are not strongly affected. 
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Figure 4. The infrared-derived star formation rate of each galaxy 
plotted against its stellar mass. The size of each plotting symbol is 
a function of the galaxy's radio spectral index, with smaller sizes 
indicating flatter spectral indices. The solid and dashed lines in- 
dicate the main-sequence and Icr (0.26 dex) dispersi on for 9 fim 
detected AKARl galaxies at z ~ as measured by lElbaz et al.l 
pOTH . SFR/M. ~ 0.25Gyr-l). Each of the infrared-luminous 
galaxies in this sample is located well above the main-sequence re- 
lation. Sources with 8.44 GHz brightness temperatures larger than 
lO**'^ K, indicating the presence of an AGN, are identified by filled 
circles. 



Like lElbaz et al.l (|20lTI ). we decide to chose the main- 
sequence relation having a slo pe of unity as this is con- 



Elbaz et al. 2007), 



(jMagdis et al.ll20iai . and 



sistent with findings at z 
dPannella et al."2009), z • 
2: ~ 4 (|Daddi et_al. 2009(1. 

We overplot our sample galaxies and find that each 
source lies more than Itr away from the main sequence, 
as expected for starbursting galaxies. Sources having the 
fiattest spectral indices tend to be furthest from the main 
sequence. 

In Figure [5] we investigate the relation between the 
distance of these galaxies from the main sequence versus 
the compactness of the starburst as measured by their 
radio spectral index. Since we have adopted a constant 
specific star formation rate to define the z ~ main se- 
quence, the distance from the main sequence is given by 
the specific star formation rate itself. We find a gen- 
eral trend where the specific star formation rate among 
the sample increases as the radio spectral index flattens. 
While a general trend does exist, the entire sample lies 
well above the z ~ main sequence. The minimum 
and mean specific star formation rates among the sam- 
ple are a factor of ~2.5 and 16 times larger than the 
z ~ main sequence value, respectively. However, many 
sources have radio spectral indices consi stent with nor - 
mal star- forming galaxies (i.e., a ~ 0.8: ICondonl[T99l . 
which comfortably lie on the main sequence. Thus, it 
seems that a fiat radio spectral index does indicate that 
a source is above the main sequence, but galaxies having 
a "normal" star-forming radio spectra (i.e., a ~ 0.8) are 
not necessarily on the main sequence. 

Since we know that typical z ~ star-forming galaxies 
lying on the main sequence have radio spectral indices 
near a ^ 0.8, it may be more appropriate to fit the data 
with something other than an ordinary least squares fit 
that passes through this normalization point. The radio 
brightness temperature, which is a measure of the surface 
brightness or "compactness" of a source, is given by the 
expression Tb = Te(l-e-^«)+Tb(0)e-^«, where Te is the 
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Figure 5. The specific star formation rate of eacli galaxy plot- 
ted as a function of radio spectral index (generally) measured be- 
tween 1.49 and 8.44 GHz (see Table [ij. Sources with 8.44 GHz 
brightness temperatures larger than 10*'^ K, indicating the pres- 
ence of an AGN, are identified by filled circles. We additionally 
identify those sources categorized as being AGN dominated by 
their low 6.2 /xm PAH EQW (circles with crosses); note, each of 
the sources with an 8.44 GHz brightness temperature larger than 
lO^-S K also has a 6.2 fim PAH EQW < 0.27 /.tm . The solid line 
is a fit to the non-AGN sources, as identified by either their low 
6.2 /im PAH EQW or high 8.44 GHz brightness temp eratu re, of the 
form sSFR = 77(1 - e°-°MS ) + sSFRms e°-°MS (see gjl. which is 
forced to fit through the typical radio spectral index (qms ^ 0.8) 
and specific star formation rate (sSFRms ~ 0.25 Gyr~^) of 2 ~ 
star- forming galaxies (diamond). The corresponding error bars in- 
dicate the Iff dispersions of these parameters (i.e., o"c«ms ~ 0.13 
and ctsSFRms = 0.26 dex) 

thermal electron temperature, Tg is the free-free optical 
depth, and Tb(0) is the initial brightness temperature. 
The specific star formation rate has been shown to be 
related to th e compactness (br ightness temperature) of 
a source (e.g.. lElbaz et al.ll2011[ ). and in Figure[5]we have 
found a relation between the radio spectral index and the 
optical depth of the 9.7 /im silicate feature. Thus, let us 
assume that the specific star formation rate is related 
to the radio spectral index with a function having the 
form of sSFR = 7/(1 - e"-""=5) ^_ sSFRms e"-""^ where 
q;ms Sind sSFRms £^re the typical radio spectral indices 
and specific star formation rates for galaxies on the z ~ 
main sequence, respectively, and 77 is a normalization 
constant related to the maximum specific star formation 
rate, which is determined by the fit to the data. 

Taking ajvis = 0-83 with a standard deviation of 
(TaMs = 0-13 (INiklas et al.lflQQTl) and sSFR - 0.25 Gyr'^ 
with a standard deviation that is a factor of 1.82 
(|Elbaz et al.l[20Tll) . we find that sources not identified 
as harboring AGN by either a low 6.2 /zm PAH EQW 
(i.e., < 0.27 //m) or high 8.44 GHz brightness tempera- 
ture [i.e., Tb(8.44GHz) > lO^-^K] are best fit by 

where the error on rj is estimated by the standard devia- 
tion of the distribution from the fit. We show this fit 
to the data in Figure [SJ along with the z ~ normaliza- 
tion point (diamond) and la dispersions. If we instead 
only exclude those sources identified as AGN by their ex- 
cessively large 8.44 GHz brightness temperatures, given 
that it is unclear whether the small 6.2 fim PAH EQW is 
completely the result of an embedded AGN, we obtain a 



slightly larger value for 77, being « 17.1 ± 4.3. 

5. CONCLUSIONS 

We have compared the mid-infrared spectral proper- 
ties of local infrared-luminous systems with high reso- 
lution radio measurements in the literature, specifically 
comparing 8.44 GHz brightness temperatures and 1.4 to 
8.44 GHz radio spectral indices with 6.2 /im polycyclic 
aromatic hydrocarbon (PAH) and 9.7 /im silicate absorp- 
tion features. In doing this, we have also investigated 
the relations between the main-sequence of star-forming 
galaxies in the local universe, as defined in the star for- 
mation rate vs. stellar mass plane, and their compact- 
ness as indicated by their radio spectral indices. Our two 
main conclusions can be summarized as follows: 

• We find that the fiattening of the radio spectrum 
increases with increasing silicate optical depth, 
8.44 GHz brightness temperature, and decreasing 
size of the radio source even after removing poten- 
tial AGN, supporting the idea that compact star- 
bursts show spectral fiattening as the result of in- 
creased free-free absorption. This trend addition- 
ally suggests that dust obscuration must largely be 
coming from the vicinity of the compact starburst 
itself and not the foreground galaxy disk. 

• We observe a trend such that the radio spectral 
index of galaxies fiattens with increasing distance 
above the main sequence (i.e. with increasing spe- 
cific star formation rate). This suggests that star- 
bursts become more compact and deeply embedded 
with increasing specific star formation rate, mov- 
ing them further away from the main sequence of 
star- forming galaxies. 

With the Karl G. Jansky VLA now online, we have 
reached a point where it is feasible to significantly im- 
prove upon the type of analysis presented in this pa- 
per by conducting higher frequency radio continuum sur- 
veys (e.g., > 10 GHz) of local st arbursts, probing <100pc 
scales (e.g.. ILerov et al.l I2011|) . For instance, by hav- 
ing resolved sizes along with better spectral coverage, to 
identify the cutoff frequency at which the free-free op- 
tical depth equals unity, one can begin to characterize 
additional physical characteristics of the starbursts such 
as their emission measures and electron densities. Ad- 
ditionally, we are likely to see a new generation of high- 
frequency, deep-field radio surveys that compliment ex- 
isting 1 .4 GHz data, while providing sub-arcsecond reso- 
lution of distant dusty starbursts. Such surveys will need 
to rely on investigations such as these to interpret their 
observations. 
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